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energy  expression,  we  discovered  a  simple  yet  efficient  de-perturbative  correction  that  eliminates  the  unphysical  perturbation 
due  to  the  finite-strength  CAP.  This  new  version  of  CAP-EOM-CCSD  has  shown  robust  performance,  improved  accuracy,  and, 
in  particular,  decreased  dependence  on  the  CAP  onset  parameter.  We  also  tested  this  approach  for  potential  energy  curves. 
Unlike  the  original  versions  of  CAP  methods,  our  1st  order  corrected  CAP-EOM-CCSD  yields  smooth  and  internally  consistent 
potential  energy  curves.  The  results  have  been  published  in  J.Chem.  Phys.  (1  paper)  and  J.  Phys.  Chem.  Lett.  (2  papers). 

We  extended  the  methodology  to  describe  properties  of  resonance  states  by  implementing  Dyson  orbitals  and  transition  dipole 
moments;  this  is  the  first  implementation  of  this  kind.  The  results  published  in  J.  Chem.  Phys..  We  also  applied  this 
methodology  to  interpret  the  experimental  results  on  resonances  in  AgF-  and  CuF-;  in  these  systems  we  discovered  new  type 
of  electronic  structure,  dipole-stabilized  shape  resonances.  The  results  were  published  in  J.  Phys.  Chem.  Lett. 

As  far  as  applications  are  concerned,  one  exciting  finding  is  that  we  have  found  a  new  stable  form  of  molecular  solid  made  of 
nitrogen  (N8).  This  might  be  of  interest  to  high-energy  compounds.  The  results  of  this  work  have  been  published  in  Nature 
Chemistry. 

Our  work  on  resonances  was  recognized  by  Phi  Kappa  Phi  faculty  Recognition  Award  in  2016.  The  nomination  to  the  award, 
which  describes  our  contributions  to  the  theory  of  metastable  states,  is  attached. 

Technology  Transfer 

The  codes  developed  under  the  ARO  support  are  integrated  into  the  Q-Chem  package.  There  were  released  in  Q-Chem  4.0. 
The  modules  are  also  available  as  open  source  (upon  request)  and  can  be  integrated  in  other  packages  (integration  in  Cfour 
and  Psi4  is  planned). 


(4)  Scientific  Progress  and  Accomplishments  (description  should  include  significant 
theoretical  or  experimental  advances) 

The  most  important  achievement  is  efficient  production-level  implementation  of  cs- 
EOM-CCSD  (EE  and  EA  variants)  and  CAP-EOM-CCSD  (EE,  EA,  and  IP).  This  allowed 
us  to  investigate  the  performance  of  these  approaches  and  to  identify  their  potential  and 
limitations.  We  have  investigated  the  performance  of  cs-EOM-CCSD  using  small  atomic 
systems;  we  have  found  that  the  strong  basis  set  dependence  of  the  results  is  due  to 
electron  correlation,  and  not  due  to  the  diffuseness  of  the  resonance  wave  functions,  as 
was  previously  thought  (this  work  was  published  in  J.  Chem.  Phys.).  Our  findings 
suggest  a  practical  way  of  developing  compact  basis  sets  for  resonances.  Specifically, 
we  intent  to  employ  Frozen  Natural  Orbitals  technique  to  develop  such  bases.  We  also 
found  that  frozen  core  method  provides  a  reasonable  way  for  dealing  with  core  electron 
problems. 

We  have  investigated  the  performance  of  CAP-EOM-CCSD  for  various  molecular  shape 
resonances.  By  analyzing  analytic  energy  expression,  we  discovered  a  simple  yet 
efficient  de-perturbative  correction  that  eliminates  the  unphysical  perturbation  due  to  the 
finite-strength  CAP.  This  new  version  of  CAP-EOM-CCSD  has  shown  robust 
performance,  improved  accuracy,  and,  in  particular,  decreased  dependence  on  the  CAP 
onset  parameter.  We  also  tested  this  approach  for  potential  energy  curves.  Unlike  the 
original  versions  of  CAP  methods,  our  1st  order  corrected  CAP-EOM-CCSD  yields 
smooth  and  internally  consistent  potential  energy  curves.  The  results  have  been 
published  in  J.Chem.  Phys.  (1  paper)  and  J.  Phys.  Chem.  Lett.  (2  papers). 

We  extended  the  methodology  to  describe  properties  of  resonance  states  by 
implementing  Dyson  orbitals  and  transition  dipole  moments;  this  is  the  first 
implementation  of  this  kind.  The  results  published  in  J.  Chem.  Phys..  We  also  applied 
this  methodology  to  interpret  the  experimental  results  on  resonances  in  AgF-  and  CuF-; 
in  these  systems  we  discovered  new  type  of  electronic  structure,  dipole-stabilized  shape 
resonances.  The  results  were  published  in  J.  Phys.  Chem.  Lett. 

As  far  as  applications  are  concerned,  one  exciting  finding  is  that  we  have  found  a  new 
stable  form  of  molecular  solid  made  of  nitrogen  (N8).  This  might  be  of  interest  to  high- 
energy  compounds.  The  results  of  this  work  have  been  published  in  Nature  Chemistry. 

Our  work  on  resonances  was  recognized  by  Phi  Kappa  Phi  faculty  Recognition  Award  in 
2016.  The  nomination  to  the  award,  which  describes  our  contributions  to  the  theory  of 
metastable  states,  is  attached. 


Metastable  electronic  states:  New  advances  on  the  old  frontier 


Figure  1:  In  high-energy  environments, 
electrons  have  sufficient  energy  to  over¬ 
come  their  attraction  to  the  nuclei  giv¬ 
ing  rise  to  the  fourth  state  of  matter 
called  plasma.  Plasma  consists  of  charged 
particles  (electrons,  ions)  and  metastable 
atomic  and  molecular  species.  If  you  think 
plasma  is  an  exotic  matter,  think  twice. 
Plasma  constitutes  99.999  %  of  the  visi¬ 
ble  universe,  as  it  is  the  main  component 
of  stars.  Plasma  is  often  generated  by 
electric  discharge  such  as  in  thunderbolt 
or  in  plasma  TVs.  Cosmic  radiation  pro¬ 
duces  plasmas  in  atmosphere  giving  rise  to 
northern  lights. 

What  are  metastable  electronic  states?  These  refer  to  atoms  and  molecules  that  cannon  hold 
all  their  electrons  and  decay  by  ejecting  an  electron,  via  so  called  auto-ionization  process.  They 
are  transient  species  that  have  finite  lifetime.  Yet,  they  play  a  crucial  role  in  science  and  technol¬ 
ogy.  These  states  are  ubiquitous  in  a  fourth  state  of  matter,  plasma,  which  can  be  generated  in 
high-energy  environment  (stars,  electrical  discharge,  cosmic  radiation,  see  Fig.  1).  They  also  appear 
when  molecules  are  exposed  to  radiation  or  interact  with  metal  surfaces,  e.g.,  in  molecular  electron¬ 
ics.  Electronic  resonances  can  serve  as  doorway  states  in  processes  such  as  catalysis,  interstellar 
synthesis,  and  DNA  damage..  They  are  also  commonly  encountered  in  experiments  using  advanced 
light  sources. 

Despite  their  importance,  theoretical  description  of  such  metastable  states  has  been  elusive.  Many 
prominent  theoreticians  have  worked  on  this  problem  since  the  dawn  of  Quantum  Chemistry.  Many 
ingenuous  ideas  have  been  put  forward,  however,  they  have  not  result  in  practical  methods  due  to 
the  complexity  of  the  many-electron  nature  of  the  problem  and  the  importance  of  electron  correlation. 

Krylov’s  group  began  to  work  on  this  problem  in  2012.  Krylov  realized  that  equation-of-motion 
coupled-cluster  framework,  a  powerful  approach  to  excited-states  and  open-shell  systems,  provides 
an  ideal  platform  for  building  the  theory  for  resonances.  The  series  of  the  5  papers1-5  present  the 
realization  of  the  approach  and  illustrate  its  capabilities.  The  essential  advance  was  to  realize  how 
to  achieve  the  stabilization  of  a  metastable  state  by  the  so-called  complex  absorbing  potential,  such 
that  its  wave  function  can  be  computed,  and  then  remove  the  perturbation  (due  to  the  complex 
absorbing  potential)  by  deperturbing  correction.3  This  lead  to  the  robust  performance  and  reliable 
results.  The  new  approach  enables  description  of  not  only  atomic,  but  also  molecular  resonances  and 
yields  potential  energy  surfaces  which  are  smooth  and  internally  consistent.4  These  tools  enabled 
the  discovery  of  the  new  type  of  resonances,  dipole-bound  shape  resonances,  which  is  presented  in 
the  last  paper  in  the  series.  The  calculations  by  Krylov  and  Jagau  revealed  that  molecules  with 
large  dipole  moments  (like  CuF  and  AgF)  can  support  shape  resonances  resembling  dipole-bound 
anionic  states.5  In  these  states,  the  extra  electron  resides  outside  of  the  molecular  core,  leading  to 
the  distinctly  different  properties.  For  example,  the  lifetime  of  such  s  shape  resonances  is  not  affected 
by  molecular  vibrations,  in  sharp  contrast  to  the  regular  shape  resonances,  as  illustrated  in  Fig.  2. 
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Figure  2:  Comparison  of  regular  shape  and  dipole-stabilized  resonances  by  means  of  potential  energy 
curves  for  the  1  2II^  resonance  of  and  the  X  ground  state  of  N2  (left,  data  from  Ref.  4)  and 
for  the  1  2II  resonance  of  AgF-  and  the  X  1X+  ground  state  of  AgF  (right,  Ref.  5).  Resonance  widths 
as  a  function  of  bond  distance  are  displayed  as  insets.  Also  shown  are  the  molecular  orbitals  occupied 
by  the  excess  electron  in  the  resonance  states.  Regular  shape  resonances  in  diatomic  molecules  such 
as  in  N2  are  stabilized  at  stretched  geometries;  the  energy  of  the  attached  state  drops  below  that 
of  the  neutral  and  its  width  becomes  zero,  which  can  be  rationalized  by  the  symmetric  shape  of  the 
molecular  orbital.  In  contrast,  energies  and  lifetimes  of  dipole-stabilized  resonances,  in  which  the 
excess  electron  resides  outside  the  molecular  core,  remain  nearly  constant. 


Individual  contributions.  In  the  5  papers  listed  below,  Prof.  Krylov  is  the  lead  theoretical 
author.  The  theoretical  work  presented  in  these  papers  has  been  conducted  at  USC  by  Krylov  and 
her  postdocs  (Bravaya,  Jagau,  Epifanovsky)  and  graduate  students  (Zuev).  The  contributions  of 
collaborators  from  Berkeley  (Head-Gordon,  Sundstrom)  and  Q-Chem,  Inc  (Shao)  are  the  computa¬ 
tional  infrastructure  (Q-Chem  base  code).  Prof.  Mabbs  and  his  group  (Dao,  Holtgrewe)  conducted 
experimental  measurements. 
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